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ABSTRACT 



Context. The two types of Fanaroff-Riley radio loud galaxies, FRI and FRII, exhibit strong jets but with different 

properties. These differences may be associated to the central engine and/or the external medium. 

Aims. The AGN classification FRI and FRII can be linked to the rate of electromagnetic Poynting flux extraction 

from the inner corona of the central engine by the jet. The collimation results from the distribution of the total 

electromagnetic energy across the jet, as compared to the corresponding distribution of the thermal and gravitational 

energies. 

Methods. We use exact solutions of the fully relativistic magnetohydrodynamical (GRMHD) equations obtained by a 
nonlinear separation of the variables to study outflows from a Schwarzschild black hole corona. 

Results. A strong correlation is found between the jet features and the energetic distribution of the plasma of the inner 
corona which may be related to the efficiency of the magnetic rotator. 

Conclusions. It is shown that observations of FRI and FRII jets may be partially constrained by our model for spine 
jets. The deceleration observed in FRI jets may be associated with a low magnetic efficiency of the central magnetic 
rotator and an important thermal confinement by the hot surrounding medium. Conversely, the strongly collimated 
and accelerated FRII outflows may be self collimated by their own magnetic field because of the high efficiency of the 
central magnetic rotator. 

Key words. MHD, Relativity, Galaxies: active, Galaxies: jets, Acceleration of particles, Black hole physics 
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1. INTRODUCTION 

According to the standard Active Galactic Nuclei (AGN) 
paradigm, their radio luminosity is related to the presence 
of powerful relativistic jets (radio loud AGN), or, to mildly 
sub- relativistic outflows (radio quiet AGN). And, by as- 
suming a supermassive BH surrounded by an accretion 
disk/torus, the different AGN phenomenologies observed 
in both classes are related to the orientation of the axis of 
the BH/disk system with respect to the line of sight, and 
the thickness of the torus which is resp onsible for the ob- 
scuration effects ( Urry fc Padovani|1995 ). Typical examples 
for radio quiet AGN are the various types of Seyfert I - II 
galaxies, where uncollimated (or loosely collimated) winds 
are outflowing from the BH/disk system at a speed of a few 
thousands km s _1 . However, besides their inclination to the 
line of sight the classification cannot be complete without 
invoking another key parameter to explain the outflow dif- 
ferences between the various AGN. These differences may 
be related to galaxy environment effects and/or intrinsic 
properties of the AGN, as shown in Fig. 1 (see e.g. Kaiser 
& Alexander 1997; Celotti 2003; Kaiser & Best 2007). 

For radio loud objects a fundamental role is played by 
Doppler boosting, strongly affecting the luminosity and 



spectral properties of these AGN. In fact, these radio 
sources are associated with powerful relativistic jets which 
reach at the parsec scale, high Lorentz factors 7 ~ 5 — 30 



(|Urry fe Padovani 1995 Kellermann et al 



2004 Piner et 



opening 



|al.||2003p . Such jets are stron gly collimatea wit 
angles of the order of ~ 3° (Pushkarev et al. 120091) . In 
nearby AGN recollimation is inferred from the inner ra- 
dio jet structure (Horiuchi et al. 2006 for Cen A; Kovalev 
et al. 2007 for M87). We recall that the main classes of 
radio loud AGN are Radio Quasars, Flat Spectrum and 
Broad Line Radio Galaxies, BL Lacs, Fanaroff Riley I (FRI) 
and Fanaroff Riley II (FRII) objects. According to the uni- 
fied model, FRI objects arc misaligned BL Lacs, while the 
parent population of FRII are Radio Quasars, Broad Line 
Radio Galaxies and the brightest BL Lacs, as sketched in 
Fig. [l] Regarding in particular the FR I and FR II di- 
choto my, we briefly outline i n the following their main prop- 
erties ( |Fanaroff fe Riley||1974| ) : 
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- In FR II sources the extended radio morphology shows a 
clear, generally one sided collimated (within a few degrees) 
thin jet, terminating into a hot spot and surrounded by 
diffuse blobs. Conversely, in FR I sources the collimated 
symmetric jets smoothly merge into the extended emitting 
regions. 
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- FR II jets look highly relativistic and narrow along their 
whole length (tens of kpc). FRI je ts are conversely relativis- 
tic only on pc scales ( Bridle|1982[), becoming sub relativistic 
and diffuse on kpc scales ( Giovannini et al.|2005 ). However, 
in some FR I sources the structure of the jet m the kilo- 
parsec scale appears more complicated, with an inner spine 
that remains relativistic and an outer shell that d ecelerates 
and becomes sub-relativistic ( Canvin et al.|[2005 |. 

- FR II sources are more powerful than the FR I ones, with 
threshold power ~ 10 25 W Hz -1 sr _1 increasing with the 
radio galaxy luminosity (Ledlow & Owen 1996). 

- FR II are usually found in poor gas environments, with 
j ets probably collimated by their helical magnetic fields 



Type 1 (Narrow Line) 



Type ] (Broad Line) 



Type (Unusual) 



(Hardcastle fe Worrall 2000 Asada et al. 2002 Zavala & 



jryl or|2005| ) and slightly interacting with the external gas 



Rich environments harbor mostly FR I sources and their 
jets, thermally confined (at least partially) and ap pearing 



to strongly interact with the intracluster medium ([Kaiser 
fe Alexander]|1997[ |Laing et al.]|1999| |Gabuzda|[2003| ) . The 



measured transverse magnetic field suggests the presence 
of intern al shocks where the tangled magnetic field is com- 



Gomez et al. |2008| ). Those 
re thermal collimation of the 



pressed ( Gabuzda et al.||1994 
shocks could be the result of t 
jet. In the following, we briefly discuss numerical versus 
analytical modeling of multicomponent jets. 

There are two main theories to interpret the above ob- 
servational characteristics. The first explains the morpho- 
logical differences as mainly due to the different physical 
properties o f the environmen t in which the relativistic jet 
propagates flDe Young|[l993[ |Bicknell||1995l |Laing et al 



1999 



Gopal-Krishna fc Wiita||2000| IMeliani et al.| [2008 ) . 
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1'he second explains the dichotomy by involving a differ 

ence in the nature of the central engine, the spin of the 
centr al black hole, the accretion rate and the jet composi- 
tion flBaum et al.||1995| [Reynolds et al.| |1996[ |Meier||1999| 
Meliani fc Keppens|2009p . Finally, there may be a combina- 
tion of external and engine factors to explain the FR I/FR 
II dichotomy, as we have suggested in Meliani et al. (2006a) 
and studied in Wold et al. (2007). 

Si milarl y to jets from Young Stellar Objects ( jFerreira 



et al.||2006[), AGN jets probably have at least two com po- 
nents ( [SoTet al.||1989[ |Tsinganos fe Bogovalov||2002[ ) one 
originating from the surrounding Keplerian disk (Baum et 
al. 1995; Meier 2002; Begelman & Celotti 2004) and the 
other from the inner corona surrounding the central black 
hole. This corona can be created by the " CEntrifug al pres- 



sure supported Bo undary Layer" model (CENBOL) (Das & 
Chakrabarti1[2002l ). The corona can also be creat ed by the 
mechanism presented in|Kazanas & Elison (1986). A third 



alternative to produce such a corona with pair plasma is the 
Blandford & Znajek (1977) model where the jet is powered 
by the spinning black hole. However observations indicate 
that the jet should have both hadronic and leptonic com- 
ponents as explained in two-component models (Henry & 
Pelletier 1991; Fabian & Rees 1995). 

The MHD equations can be solved through numerical 
simulations, which describe the evolution of the jet config- 
uration. The availability of more and more powerful com- 
puting facilities and sophisticated numerical codes allows 
a quite complete descrip tion and understanding of the jet 
acceleration/collimation ( |Komissarov et aL"[|2007| [ Porth fc 
Fcndt 2010) and the accretion/ejection process (Koid 



ndt||2010[) and the accretion/e j ection process QKoide et 
||1998[ |1999| |McKinney||2006l |McKinney fc Blandford 



2009 Gracia et al. 2006[ |2009[ ). Recent numerical Simula- 
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Fig. 1. Standard classification of AGN sources following 
Urry & Padovani (1995). The horizontal axis represents the 
inclination of the source axis with the line of sight. The ver- 
tical axis we suggest that it may be linked to the efficiency 
of the underlying magnetic rotator to collimate the flow. 



tions have progressed to general relativis tic magnetohydro- 
dyn amic (GRMHD) jet launching, as in McKinney (2006) 
and Hardee et al. (2007), suggesting also the formation of 
jets with two components. However computational limits 
do not allow yet to follow simulations for very long times 
and reach exact stationary configurations. It also fails at 
analyzing structures with very different scale lengths. 



Nevertheless, tremendous progress on understanding 
the physics of relativistic jet acceleration/deceleration - 
and therefore the FRI/FRII dichotomy - has been done 
thanks to numerical simulations of jet propagation in 
the asymptotic regions. Some investigated the relativis- 



lar medium ( 


Duncan & Hughes 


1994 


Marti et al. 


1997 


Komissarov & Falle||1998| |Aloy et al. 


1999| |Rossi et al. 



FR II jets may be a consequence of the power of the jet. 
Many groups had also investigated the two-dimensional rel- 
ativistic magnetized jet propagation in an external medium 



van Putten||19 96| 
Keppens et al 



6) Komissarov|[l999 Leismann et al.|[2005 
2008P and 3D (e.g. Mizuno et al. 2007, 



Mignone et al. 2009). They showed that the interaction 
between the jet and the external medium depend on the 
magnetization of the jet and the density ratio between 
the jet and the external medium. The role of the envi- 
ronment may be crucial in HYbrid MOrphology Radio 
Sources (HYMORS), as shown by Gopal-Krishna & Wiita 
(2000, 2002). These radio sources appear to have a FR II 
type on one side and a FR I type diffuse radio lobe on 
the other side of the active nucleus. This last model for 
HYMORS has been recently confirmed by numer ical simu- 
lations of two component jets ( |Meliani et al. 2008). Another 
alternative to the FRI/FRII dichotomy could come from 
the nature of the instabilities that develop in the jet, or, 
at the jet interface with the external medium, or, with 
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anoth er surrounding outflow component (Keppens et al. 
2008). This is consistent with the fact that hydrod ynarrT 
ical jets wi th high Lorentz factors are more stable (Marti 
et al. 1997|). Poloidal magn etic fields also help to stabilize 
the jet ( |keppens et al. 2008[). On the last alt ernative note 
that simulations by Meliani "Sj Keppens (2009) confirm that 
the deceleration in FR1 jets could be attributed to a strong 
Rayleigh Taylor instability between the spine jet and the 
surrounding component, supporting a two component jet 
structure. 

In parallel to the development of time dependent simu- 
lations, steady jet solutions in GRMHD were first obtained 
numerically by solving the transfield equation in the force- 
free limit (Camenzind 1986). This study has been further 
developed in GRMHD by using first a Schwarzschild met- 
ric and then extending it to a Kerr metric (Fendt 1997). 
However this method cannot incorporate consistently the 
mass loading of the jet. 

Exact models for disk winds can be constructed via a 
nonlinear separation of the governing MHD equations. This 
technique of radial self-similarity has been largely explored 



thors ( 


Bardeen & Berger 


1978 


Blandford & Payne 


1982 


Li et a 


1992 Contopoulos 1994 


Ferreira 1997; Vlahakis & 


Konigl 


2004). 



aration which is used to produce models of pressure driven 



winds ( Meliani et aTj 2006a I . It is also complementary to 
magnetically driven disk winds. Such models may describe 
the inner spine jet from the central object necessary to sus- 
tain the outer disk wind but where the radial self-similar 
models fail by construction. In this paper, using such solu- 
tions, we show that the collimation criterion developed in 
the frame of this model can help understanding how the 
FRI/FRII dichotomy may influence the morphology of the 
inner spine jets. As a back reaction, the spine jet dynamics 



2009). 



Fendt 


2009(), even if 


eak m 


^atsakos et al. 



In Sect. [2] we recall briefly the main assumptions of the 
model. In Sect. [3] we summarize the details of the standard 
AGN jet classification and how it helps to constrain the 
parameters. In Sect. [4] we present an interesting solution 
for FRI type spine jet and in Sect. [5] another one for jets 
associated with FRII objects. In Sect. 6 we discuss and 
summarize the main implications of our model. 

2. Model assumptions 

2.1. The self-similar model 

We use the ideal GRMHD equations in the background 
spacetime of a Schwarzschild black hole, and neglect the 
effects of self gravity of matter outside the black hole. The 
spacetime curvature at a distance r from the black hole is 
given by the lapse function, 



h = a 1 



r 



(1) 



where vq is the Schwarzschild radius. 

Following Meliani et al. (2006a), all physical quantities 
are normalized at the Alfven radius r* along the polar axis 
where the meridional angle is zero (9 — 0). We define a 



dimcnsionless spherical radius R — r/r+ 
G, and magnetic flux function a, 

R 2 . 9 „ „ r sin(6») 



G 2 {R) 



sin 



G = 



r* sin(0*) 



cylindrical radius 



(2) 



To describe the GRMHD outflow of the coronal plasma, 
we use the relativistic meridionally self-similar solutions 
presented in Meliani et al. (2006a). The specific enthalpy, 
and density in the lab frame, together with the pressure, 
velocity and magnetic field are given in terms of functions 
of the radial distance R, 



1 



h~fw = 7j*7*w* ( 1 - 
P 

V r 
Vn 



Po + -7*V^rV^n(i?)(l + net) , 



V+M 2 



1 



h 2 G 2 y/i + 5a 
V+M 2 hF 1 



M 2 7Vb 
v 2 D 



Va, = - 



h 2 G 2 2 VTTfo 
h AK N v Rsm9 



sin 9 . 



K G 2 D VTTfo 



G 2 



cos 9 . 



B r = 
B e = 
B v = 
where 



N B =*-G 2 , 



B*hF 



\B+ K N B 
G 2 h D 



i?sin( 



M 2 h 2 



M 2 



(3) 
(4) 
(5) 
(6) 

(7) 

(8) 

(9) 
(10) 

(11) 
(12) 



The free parameters 5 and k describe the deviation from 
spherical symmetry of the ratio of number density /enthalpy 
and pressure, respectively, while A is a constant controlling 
the angular momentum extracted by the jet. The constants 
v and fi measure the escape speed in units of the light speed 
and the escape speed in units of Alfven speed, at the Alfven 
point along the polar axis, respectively, 



V'i 



V m 



fi 



(13) 



Thus, all physical quantities arc determined in terms of 
constant parameters, S, k, A, n, v and the three unknown 
functions, n(i?), F(R) and M 2 (R). Note that Il(R) is the 
dimensionless pressure function, defined modulus a con- 
stant Pq, F is the expansion factor, and M is the poloidal 
Alfvcnic number, 



F 



dlnG 2 
dlni? 



M 



Airh 2 nw^ 2 V 2 



B 2 c 2 



(14) 



These three unknown functions Tl(R), F(R) and M 2 (R) are 
determined by three nonlinear equations. We start integrat- 
ing these equations from the Alfven critical surface, tak- 
ing into account there the corresponding regularity condi- 
tion and then integrate downwind and upwind, crossing the 
other critical points, for details see Meliani et al. (2006a). 
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The light cylinder is defined by the function x = CIL/E 
becoming unity, where f2 is the angular speed, L the total 
angular momentum per unit mass and E the generalized 
Bernoulli integral. It is a measure of the energy flux of the 
magnetic rotator in units of the total energy flux. The s-s 
description is possible only if the jet is rotating at subrela- 
tivistic speeds. In such conditions x must remain small and 
therefore the light cylinder effect have to be negligible for 
our solution to be valid. 

In this setup of the relativistic MHD problem an extra 
parameters exists e which is constant everywhere (Meliani 
eTaT|2006a| ): 



hfR 2 G 2 



F 2 

T 



i 

V 2 



R 2 
'"h 2 G 2 



G 2 h 2 



D ) 



{5 - k) v 2 
h 2 R 

2A 2 N B 



(15) 



This is the relativistic generalization for a Schwarzschild 
black hole of the classical constant found in Sauty et al. 
(2004) that measures the magnetic energy excess or deficit 
on a nonpolar streamline, compared to the polar one. To 
first order, e determines the fraction of Poynting flux car- 
ried by the jet in the asymptotic region, such that 1 — e 
measures the fraction of the Poynting flux which is used 
in the jet acceleration. Thus, if e > we have an Efficient 
Magnetic Rotator (EMR) where magnetic collimation may 
dominate, while if e < we have an Inefficient Magnetic 
Rotator (IMR) where collimation cannot be of magnetic 
origin but should be thermal, if the flow is collimated. 

3. New interpretation of the Fanaroff-Riley 
classification 

We propose here to examine the vertical classification in 
Fig. [l]by means of our parameter e, or less ambitiously, to 
interpret how the various observations we have on FRI and 
FRII jets may influence the formation of the inner spine 
jet co mponent. As in the non relativistic case (Sauty et al. 
1999), this parameter e allows to classify jets according to 
the efficiency of the central magnetic rotator. Jets emerg- 
ing from efficient central magnetic rotators, e > 0, colli- 
mate cylindrically without oscillations in the asymptotic 
region. In this type of jets, the velocity increases monotoni- 
cally to reach its asymptotical maximum value. Conversely, 
jets associated with inefficient central magnetic rotators, 
e < 0, are collimated mainly by the pressure of the exter- 
nal medium. Thus, this type of jets strongly interacts with 
the ambient medium and this induces oscillations in their 
shape at the asymptotic region. The speed of these jets also 
does not increase monotonically but it oscillates too. The 
plasma in those jets is accelerated until an intermediate re- 
gion where the speed reaches its maximum value. Then, the 
acceleration of the jet stops in this region. Further away in 
the recollimation region the outflow slows down. Note that 
such jet oscillations could lead finally to a more turbulent 
outflow consistently with the numerical simulations of two 
component jets we mentioned in introduction. After the rec- 
ollimating region we cannot exclude the presence of shocks 
or instabilities. Thus it is difficult to know if the oscillations 
that appeared in the solutions would be observable or not. 

In summary, we propose that a difference between these 
two types of spine jets associated with FR I and FR II radio 



galaxies may result from the competition between the mag- 
netic and thermal confining mechanisms. Magnetic pinch- 
ing and pressure gradient tend to compensate the trans- 
verse expansion of the jet because of the centrifugal force 
and the charge separation which also induces an outwards 
electric force. These two expanding forces are characterized 
by the free parameter A 2 , that measures the quantity of an- 
gular momentum carried along the streamlines. Magnetic 
collimation is controlled by the parameter e/A 2 (see Sauty 
et al. 1999). On the other hand, thermal collimation is con- 
trolled by the parameter n that defines the transverse vari- 
ation of the pressure. 

Hence, an appropriate manner to classify the different 
jet solutions according to the nature of their collima tion is 
based on the two free p arameters k/2A 2 and e/2A 2 (Sauty 
et al.||1999[ |2002[ [2004}. The higher is the value of those 



parameters, the stronger is the collimation and the lower is 
their terminal speed. This can be seen in Tab. [3] where we 
plot for various solutions the asymptotic jet speed in units 
of the escape speed from the base of the corona: as the effi- 
ciency of collimation increases, the efficiency of acceleration 
decreases. In other words, tightly collimated jets (larger val- 
ues of e/2A 2 , or, k/2A 2 ) have lower terminal speeds. Thus 
the strong interaction with the external medium would nat- 
urally result into a decelerated flow even if it remains stable. 
Parallely, the jet being denser would radiate more on large 
scale before the terminal shock with the ambient medium. 
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Fig. 2. In (a) we plot the Alfvcn number M 2 and, in 
(b), the cylindrical jet cross section G 2 as functions of 
the distance for k/2X 2 = 0.05 and four values of e/2A 2 
(-0.5,-0.5 x 10" 2 ,0.5 x 10~ 2 ,0.5), corresponding to Tab. 
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1.07 
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30 


13 
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90 


75 


65 


60 


10 


1.02 



Table 1. The asymptotic jet speed in units of the escape speed from the base of the corona, for various values of the 
parameters k/2X 2 and e/2A 2 



Tab. [3] also shows that jets from inefficient magnetic ro- 
tators are more powerful in transforming thermal energy 
into kinetic energy, than those from efficient magnetic ro- 
tators. There are two reasons for this. 

First, in a EMR flow, the centrifugal force at the base 
of the jet is important. Then, the last stable orbit of the 
plasma gets closer to the central black hole. Therefore, the 
corona extends closer to the black hole horizon. Thus, if 
the available total amount of thermal energy is the same 
in the corona, the plasma gravitational potential increases. 
Consequently, in EMR as e/(2A 2 ) increases, more thermal 
energy is tapped in order to allow for the plasma to escape 
and less is left for accelerating it. 

Second (see Fig. [2 1, an increase of the magnetic rota- 
tor efficiency limits the initial expansion of the outflow. 
The pinching magnetic force gets stronger after the Alfvcn 
surface and the conversion of thermal energy into kinetic 
energy stops when the jet reaches its asymptotic cylindrical 
shape. This decrease of the jet asymptotic speed with the 
increase of the magnetic rotator efficiency may seem contra- 
dictory with the usual picture. However, in axial outflows 
the contribution of the Poynting flux to the total accelera- 
tion remains weak. This is of course different from relativis- 
tic disk wind models where the acceleration is dominated 
by the conversion of Poynting flux to kinetic energy flux 



(iLi et al|1992"l |Contopoujos|1994[ |Vlahakis fc Konigl|2 004 ) 

Those jets are characterized by a strong inclination of the 
magnetic field lines and the rotation at the base of the out- 
flow is almost Keplerian. Note also that, as expected, the 
efficiency of the acceleration also increases with the degree 
of expansion, i.e., as k/2X 2 decreases. 



3.1. Estimation of the free parameters 

In the following sections, we mainly explore two examples of 
typical solutions for relativistic jets. One is associated with 
an IMR such that the contribution of the external pres- 
sure is comparable to the magnetic one, while the second 
is associated to an EMR such that the jet is self collimated 
magnetically. In order to find these solutions we first make 
an estimation of the free parameters of the model using our 
knowledge of the properties of FRI and FRII jets, in par- 
ticular in the launching region on one hand and in the far 
asymptotic region on the other hand. 

At the base of the outflow, we assume that the corona 
starts at the radial distance of the last stable orbit which is 
around the radius of the Schwarzschild black hole, as dis- 
cussed below. We consider that the last open streamlines, 
which emerge from the corona, should be in sub-Keplerian 
rotation if they are anchored in the thick disk surrounding 
the central black hole. The last open streamline in the coro- 
nal jet is the one that crosses the equatorial plane at the 
edge of the magnetic dead zone, which we assume to have 
a dipolar configuration. 



We use also some observational constraints in the 
asymptotic region of the jet. The opening angle of jet has 
to be a few degrees. Six degrees is the value inferred for the 
well measured jet of M87 (Biretta et al. 2002; Kovalev et 
al. 2007). We guess that the spine jet opening angle is even 
smaller and we took a value of ~ 1° at 1 pc. Then, from the 
expression of the dimensionless magnetic flux a, by know- 
ing the last streamline, we can deduce the asymptotic value 
of Goo ■ 



Goo — 



: sin f 



(16) 



where a cx t is the last open streamline in the jet. As a 
reasonable estimate, we have taken a cx t = 4. In fact, for 
smaller values of a cx t, the value of Goo becomes too large 
and therefore unrealistic for our model. We further assume 
that in the asymptotic region the Alfven number is of the 
order of Moo — 5: this is the order of magnitude found in 
the literature for relativistic magnetohydrodynamical jet 
propagation (Leismann et al. 2005, Keppens et al. 2008). 

We know that the asym ptotic Lorentz fac tors should 
be between ~ 3 and 10 (see Piner et al. 20031. From the 
expression of the asymptotic velocity along the polar axis 
(Eq. [6]) we can deduce from Eq. 13 the value of the free 
parameter v 2 ', 



(17) 



We used on purpose a rather lower limit for the Lorentz 
factor in order to avoid large effective temperatures in our 
model. In fact, we can use the same solutions and scale them 
up to obtain higher Lorentz factors but then the effective 
temperature would attain extremely large values above the 
mass temperature. However the corresponding high pres- 
sure could have a large contribution from a turbulent mag- 
netic or ram pressure component in the jet (see Aibo et al. 
2007, for the solar wind). In such a case the kinetic temper- 
ature would be lower. Nevertheless, as far as the collimation 
is concerned this does not affect qualitatively our discussion 
on the dichotomy between FRI and FRII and we kept this 
relatively low Lorentz factor. 

We used also the observed mass loss rate in the outflow 
to constrain the free parameter S. However, the spine jet 
probably carries only a small fraction of the observed en- 
ergy flux in AGN jets, Lj ct ,Kin ~ 10 43 ergs/s (Allen et al. 
2006) The energy flux of the coronal wind remains weak 
compared to the total mass carried by the disk-wind which 
is supposedly denser ( Vlahakis & Konigl 2004) , a situation 
similar to stellar jets associated with Young Stellar Objects 
(| Meliani et al.|2006b[ ) . Then S is deduced from the assumed 
value for the mass loss rate, 



M = 2 



^part h 7 n V p dS , 



(i 
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= tt 2 -Vll n *w* c / \ — r^da. (18 

M ' ^ Jo V w ' 

We get an equation for <5, 

(1 + a cxt 5) 4 x C* 2 - (2C 2 + 1) x (1 + a cxt <5) 3 + 
{l + a cxt 6) 2 x Cx (2 + C) 

-2xCx (l + Oext*) + l = 0, 

where C is another constant given by, 



(19) 



C = 



3/x 3 /V(l - /Lt)M 



87rr|mp art \/l - r G /r 00 7 00 n 00 M ! ^ao 



(20) 



The variable n^ is the asymptotic density and m par the 
average mass of the particles. We consider a proton-electron 



fluid, m 



par 



'"'proton ■ 



For the mass loss rate, we choose M = no 10 6 A/Edd, 
where MecM is the Eddington mass loss rate. It corr esponds 
to value fou nd for the relativistic Parker wind (Meliani 
et al. 2004|. The asym ptotic density is tak en equal to 



x no cm 3 ( jMeliani et al. 2004), with no a 



dimcnsionlcss free parameter. 

We suppose that the rotation is sub-Keplerian on the 
last open streamline of the jet at the equator. The param- 
eter r\ measures the deviation of the rotation function SI 
from its Keplerian value, 



SI = (1 - r?) 



2rl 



o/< 



(21) 



From the definition of SI (Eq. 21 ) wc can deduce the value 
of the free parameter A which is the constant controlling 
the angular momentum extracted by the jet, 



A = (i — n) 



2^ h^Jf 



(22) 



As we mention earlier we assume that the corona forms 
above the last stable orbit at tq — 3ra- We choose a typi- 
cal magnetic lever arm (i.e. Alfven radius) of 10 times the 
Schwarzschild radius, r± = lOrg. It gives, 



Ro = — = 0.3 . 



(23) 



The parameter k, which is the relative variation of the pres- 
sure with latitude, is calculated from an approximate ex- 
pression of the coronal base where M — > and dM^ /dR is 
finite in Eq. A. 2 - A. 5 (Meliani et al. 2006a) using, 



2A 2 

k = 6- Rq— (1 - /It) 



(24) 



To summarize, our free parameter estimates are: 
1. The Alfven radius 

r* = 10 xr G ^/» = 0.1. (25) 



2. The jet opening angle at ~ 1 pc is ~ 1° which 
gives the asymptotic value of Goo for the last streamline 
«cxt = 4. 

3. For the asymptotic Alfven number we chose ~ 5. 

4. The asymptotic Lorentz factor is taken to be 7oo ~ 2. 

5. The asymptotic density is taken equal to rirx, = uq x 
10~ 8 cm" 3 . 

6. The new parameter r\ which measures how sub- 
Keplerian is the velocity (the values differ from solution 
to solution). 

7. The corona is supposed to be formed above the last 
stable orbit, r c = 3rc which gives k. 

The parameter r/ is not independent from the rest of 
the model. However it measures precisely the rotation of 
the footpoints and controls the efficiency of the magnetic 
rotator to collimate the jet. Thus we have an indirect way 
to determinate the properties of the disk from the asymp- 
totic characteristics of the jet. If 77 — !► open streamlines an- 
chored in the accretion disk rotate with an almost Keplerian 
profile. On the other hand, if n — ¥ 1 open streamlines an- 
chored in the accretion disk rotate very slowly and the 
Poynting flux injected in the jet is low. 

Before going further into our modeling of FRI and FRII 
jets, we need to introduce the notion of effective tempera- 
ture. 



3.2. Temperature 



In polytropic relativistic winds (Meliani et al. 2004), the 
temperature is usually defined by the ideal gas equation 
of state P/n = ksT. Therefore, the knowledge of specific 
pressure P/n, the specific thermal energy e t h and the den- 
sity describe completely the thermodynamics of the fluid 
(temperature T, enthalpy w and pressure P). However, as 
discussed in Meliani et al. (2004), flows cannot be adia- 
batic. The polytropic approximation is just a convenient 
way to mimic heating going on in the flow. Therefore pres- 
sure, temperature and enthalpy are not the real ones but 
effective quantities that hide the extra necessary heating. 

The temperature definition in meridional self-similar 
models is similarly delicate. In fact, as indicated in Sauty 
& Tsinganos (1999) and Meliani et al. (2006a) the total gas 
pressure is not necessarily limited to be the kinetic pressure. 

Moreover, the generalized specific thermal energy of the 

p 

model, eth = (w — mc 2 ) , is not restricted to the thcr- 

v ' n 

mal energy. In self-similar model these two quantities are 

also effective quantities. They account for different physical 
processes of energy and momentum transport and dissipa- 
tion. They can include the contribution of magnetohydro- 
dynamic waves and viscous and/or radiative mechanisms. 
As a matter of fact, the complexity and variety of MHD 
processes that can contribute to the internal energy of a 
magnetized fluid, makes the definition of the real thermal 
energy impossible. 

Therefore the quantities T c g = P/n and e t h.cff = eth are 
simply the specific effective temperature and thermal en- 
ergy imposed by the dynamics of the outflow. They do not 
necessarily represent neither the kinetic temperature nor 
the thermal energy. However they are simple tools to ana- 
lyze the energetics of the flow. In the following, we discuss 
the thermodynamical properties of the fluid with this ef- 
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Fig. 3. Plot of the morphology of the solution correspond- 
ing to a FRI-type spine jet. In (a) is shown the projec- 
tion of the streamlines on the poloidal plane. The solid 
lines in the center correspond to lines where the conditions 
x\G 2 < 1CT 2 and (2 + 5a) /{I + 6a) 2 -2 < 1CT 2 are sat- 
isfied. The dashed lines correspond to x\G 2 < lO^ 1 and 
(2 + 5a)/ [1 + 5a) 2 -2 < 1CT 1 . The dashed-dotted lines cor- 
respond to x\G 2 > 1CT 1 and (2 + 5a) /(l + 5af-2 > KT 1 
(see Meliani et al. 2006a for details). 



fective temperature and not the kinetic temperature which 
we cannot calculate. 



T, 



. a = t — = ^7-7* w*m n 



(26) 



eth.eff 



w — mc 



P 

n 

1 Q W+ O 



aa 
5a ' 



(27) 



Thus along the polar axis, the specific thermal energy 
is deduced from the Bernoulli equation as follows, 



eth.eff = w*- 



V* 7 
ho 7* 



mc 



i 7 * 2 «a/ 2 n, 



(28) 



We also define Q the heat content added to the fluid, which 
is the difference between the total effective internal energy 
of the fluid eth,cff and the internal energy of the fluid ob- 
tained if it were adiabatic: 



eth.eff 



UlC 



K pol ri3- 



(29) 



where K po i is a constant. It is determined from pressure 
and density at the flow boundary either in the asymptotic 
region or in the launching region (in adiabatic flows, we 
have Kpd = P/ns). 

4. Application I - Model of FRI spine jets 

4.1. Parameters 

The first solution we show here is adapted to model the 
spine jet of radio-loud galaxies of FRI type. The environ- 
ment of such jets, i.e. the host galaxy, is known to be rich 
and containing dense gas. Moreover, the properties of FRI 
jets on the pc scale are quite different from that on kpc 
scales. In fact, in the region close to the nucleus (on the 
scale of a pc), FRI jets are accelerated to highly relativistic 
speeds. 

Beyond this region, the jets interact with the external 
medium which is denser. This interaction induces an ob- 
served deceleration of the jet. Thus, we assume that the 
outflow is likely confined by the pressure of the ambient 
medium, at least partially. In our model, these types of jets 
correspond to solutions associated with inefficient central 
magnetic rotators. In this solution 77 = 0.90 which gives an 
inefficient magnetic rotator. 

The solution corresponds to the following parameters, 



(30) 



4.2. Morphology of the FRI-type jet 

As seen in Fig. [3j the jet solution shows an initial expan- 
sion up to a distance of 100 Schwarzschild radii which then 
stops and the jet recollimates. The expansion of the jet is 
due to the strong initial inertia of the plasma carried along 
the external streamlines. The jet becomes collimatcd once 
it interacts with the "external" ambient medium that com- 
presses it. What we call "external" medium refers to the 
gas surrounding the last valid streamline of the solution; 
this can be the actual external medium of the host galaxy, 
but considering the transverse expansion of the solution, 



/' 


= 0.1, 




= 0.69, 


A 


= 0.85, 


6 


= 1.45, 


K 


= 0.42, 


e 


= -0.8475 
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it more likely corresponds to the over-pressured gas of the 
surrounding disk wind. 

The jet compression for inefficient magnetic rotators 
generates strong oscillations in the asymptotic region o f 



the jet even in the relativistic case ( Meliani et al 
The light cylinder of this jet solution is at infinity. 



2006a) 
These 



oscillations results from a transfer of energy between the 
enthalpy and the Lorentz factor as h'jw remains constant 
to first order with respect to a. 



J 




10 



r/r_ 




Fig. 4. In (a) we plot the energetic fluxes normalized to 
the mass energy of the first solution for a " FRI-type" spine 
jet. In (b) we plot for the same solution, the Lorentz factor 
7 along four different streamlines. The solid line correspond 
to the polar axis and the dotted line to the last streamline 
connected to the central corona. Other lines are intermedi- 
ate ones. 



4.3. FRI-type Jet kinematics 

The jet acceleration occurs mainly in the intermediate re- 
gion (Fig. |4b) where gravity becomes weak. Thus, all the 
enthalpy still remaining in the outflow is converted to ki- 
netic energy. 

The maximum Lorentz factor in this solution (7 = 2.8) 
is obtained at the largest expansion of the jet radius, 
where the effective temperatures are already ultra relativis- 
tic (« 10 13 K). To explain such a high effective tempera- 
ture, part of the pressure must be of non kinetic origin, with 
contribution from turbulent ram and magnetic pressure, as 
we discussed in Sec. 3.1 Thus this value of 7 is close to the 
lower observed ones, but already corresponds to a highly 
turbulent medium. The model can produce higher Lorentz 
factors provided the turbulent pressure level is sufficient. 
This is not the main topic of the qualitative discussion we 
address on the collimation of the jet itself. 



In the region of re-collimation, the increase of the pres- 
sure induces a deceleration of the jet. The Lorentz factor 
decreases from 7 = 2.8, its maximum value before the re- 
collimation, to 7 — 2 in the asymptotic region. As a matter 
of fact such a deceleration is indeed a characteristic of FRI 
jets as we mentioned. It is remarkable that this solution 
shows clearly that the main effect of the re-collimation by 
some external pressure is a global deceleration of the out- 
flow, as observed in FRI jets. The distance of recollimation 
in this solution is however smaller than the usual parsec 
scale. This may be due to the fact that here we are dealing 
only with the inner part of the jet, while observations may 
correspond to the surrounding disk-wind, or, it can be due 
to the fact that our Lorentz factor is too low. Moreover, we 
assume here that the external pressure of the host galaxy 
somehow is transmitted to the disk-wind, which in turn 
confines the spine jet. Keeping in mind these necessary pre- 
cautions, this result by itself seems interesting, if we take 
into account the simplicity of our model. 

4.4. FRI-type Jet energetics 

The temperature profile is characterized by four different 
regimes. The three first are common with the "FRII" jet so- 
lution and we shall discuss it later on. The fourth one corre- 
sponds to the asymptotic recollimated region of the "FRI" 
jet. There, the effective temperature reaches the high value 
of T e g ~ 10 13 K, because of the strong compression of the 
outflow by the external medium (Fig. [BJd). This effective 
temperature is high compared to the observed temperature 
in AGN jets that are usually of the order of T - 10 8 K. This 
large difference can be explained however by some increase 
of the contribution of non thermal mechanisms to the effec- 
tive specific thermal energy eth,off- As we mention before, 
we cannot directly compare this effective temperature with 
the kinetic temperature in the frame of this model. 



5. Application II - Model for FRII spine jets 

5.1. Parameters 

Conversely to the case of FRI, the environment of FRII 
radio-loud galaxies is relatively poor. Thus, jets from FRII 
galaxies should interact only slightly with the ambient 
medium. In these outflows, the velocity increases contin- 
uously until the asymptotic region. In fact, unlike FRI out- 
flows, the velocity in FRII jets is relativistic both on the 
parsec and kilo-parsec scales. Besides, FRII jets are so well 
collimated on large scales that they are very likely to have 
an asymptotic cylindrical shape. In our model these types 
of jets correspond to a solution associated with an efficient 
central magnetic rotator. Therefore, the value of the free pa- 
rameter rj is probably smaller, i.e. the rotation is closer to 
Keplerian velocity than in the case of the solution for FRI. 
This increases the available Poynting flux at the base of 
the jet. Therefore, the jet will be collimated by the toroidal 
magnetic pinching without any oscillations. 
The parameters for this specific solution are: 



/' 

A 
6 



0.1, 
0.65, 
1.051, 
1.35, 



(31) 
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Fig. 5. In (a) we plot the temperature profile and in (b) 
the density profile for the FRI-type jet solution. Density 
is normalized to no such that the mass loss rate is M = 



k = 0.3, 
e = 0.334. 

This is an Efficient Magnetic Rotator. It corresponds to 
a slightly lower value of r] = 0.86. This shows that the 
efficiency is very sensitive to the variations of the rotation 
frequency. 

This FRII-type jet solution is characterized by a contin- 
uous expansion up to a distance of about 100 Alfven radii 
again, but the outflow after this distance remains cylindri- 
cal, slightly expanding further out. This expansion is re- 
lated to the high magnetic pressure at the base, together 
with the strong gravity, in addition to some non negligible 
contribution of the force of charge separation. However, the 
collimation of the jet in the asymptotic region is exclusively 
of magnetic origin. It is induced by the toroidal magnetic 
pinching force and the transverse magnetic pressure. These 
two forces balance the centrifugal and charge separation 
forces. Oscillations in the jet are very weak, due to the rel- 
atively small contribution of the thermal confinement, as 
expected. 

The opening angle of the last open streamline of the 
solution at a distance of one parsec is only 0.1°, which is 
rather small compared to our initial guess. This definitely 
rules out the possibility to describe the whole jet uniquely 
with this model. Instead, we prefer to see it as the spine 
or inner part of the jet that carries away the angular mo- 
mentum of the central black hole. The situation is similar 
to what happens in Young Stellar Objects where the stel- 
lar jet is responsible for the spinning down of the protostar 
while the outer disk wind is responsible for the observed 
mass loss. 



(a) 



Light Cylinder 




Equatorial axis 



(In 



Streamline 




Magnetic Held 



Fig. 6. The same as in Fig. 3 for the second solution cor- 
responding to a FRII-type spine jet. 



5.2. FRII-type jet kinematics 

The acceleration of this solution is continuous. First, there 
is a small but effective thermal acceleration in the lower re- 
gion of the corona. In this region, the high thermal energy 
both accelerates the fluid up to 0.4c at a distance of 6rc and 
enables it to escape from the deep gravitational potential. 
A second stronger thermal acceleration of the jet occurs 
beyond the Alfven surface up to the collimation region. In 
this region, the pressure drops rapidly and asymptotically 
goes to negligible values. Therefore, the thermal energy is 
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<x = 0.0249 
cc = 0.012 



ot = 0.037 




Fig. 7. In (a) we plot the energetic fluxes normalized to 
the mass energy of the second solution for a "FRII-type" 
spine jet. In (b) we plot for the same solution, the Lorentz 
factor 7 along four different streamlines. The solid line cor- 
respond to the polar axis and the dotted line to the last 
streamline connected to the central corona. Other lines are 
intermediate ones. 



transformed into kinetic energy more effectively. In this in- 
termediate regime the velocity in the flow increases from 
0.4c to 0.92c on a scale of the order of 200rc. 

5.3. The light cylinder of the FRII-type solution 

The light cylinder in the asymptotic region of the jet (Fig. 
[6| is roughly vertical and asymptotically parallel to the 
poloidal streamlines, that remain inside the light cylinder. 
However, in this solution conversely to the previous one, 
the light cylinder is not at infinity but at a distance of 
about 20 tq from the polar axis. This reduces the domain 
of validity of the solution around its axis (cf. Meliani et al. 
2006a) where the effects of the light cylinder can be ne- 
glected. In fact self similar disk wind models can produce 
solutions crossing the light cylinder (Vlahakis & Konigl, 
2003a, b). Such solutions undergo a strong magnetic accel- 
eration ideal to obtain high Lorentz factor in GRBs for 
instance. This is of course not necessarily the case for the 
spine jet which can be accelerated by other means than the 
magnetocentrifugal process. 

5.4. FRII-type jet energetics 

The effective temperature profile of this second solution 
goes through four different regimes. 

First, in the lower corona, the effective temperature in- 
creases extremely rapidly (Fig. pa) from about 10 12 K at the 
base up to about 3 x 10 12 K. This increase is due to some 
strong initial heating in the expanding corona. The large 



Fig. 8. In (a) we plot the temperature profile and in (b) 
the density profile for the FRII-type jet solution. On the 
left the vertical lines delimitate the various domains of the 
temperature profile. On the right, density is normalized to 
no such that the mass loss rate is M = nolO -6 A/Edd- 



expansion induces a strong decrease of the density, but be- 
cause of the heating, the pressure decreases less rapidly. 

Second, in the intermediate region, the effective tem- 
perature still increases up to its maximum value of about 
3 x 10 12 K, after a relatively small decrease, because of the 
global expansion and drop of the density. 

Third, we have a transition region after the maximum 
and the asymptotic part. The effective temperature de- 
creases again to attain values around 10 12 K. This decrease 
is induced by the magnetic compression of the jet that 
brakes the density decreasing (Fig. [8}d). 

The effective temperature obtained in this solution, is 
of course high compared to observed temperatures in AGN 
jets which is typically of the order of 10 8 K in the asymptotic 
regions. As in the case of FR I this indicates that the 
contribution of non thermal energies to the acceleration and 
the heat of the jet play a relevant role. 



6. Conclusions 



We h a ve app lied exact GRMHD solutions from Meliani 
|et al.| ([2006a ) to the can onical classification of AGN jets 
( |Urry & Padovam 1995) according to their morphology. 
Our model is constructed in the frame of general relativity 
using the metric of a central Schwarzschild black hole. It 
validates the classification of AGN proposed in Sauty et al. 
(2001) where the classical MHD solutions were used with 
some additional features due to relativistic effects. In this 
study it was proposed that the inner regions of jets (spines) 
are collimated by an external denser medium in FRI and 
by the force of magnetic pinching in FRII. 
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We first proposed a method to estimate the free pa- 
rameters of the model from the known properties of AGN 
jets. In particular, the departure from Keplerian rotation 
at the footpoints of the fieldlines is measured by an extra 
parameter 77. It turns out that we get very different classes 
of solutions by slightly changing this parameter. 

First for 77 = 0.90 we obtained a recollimating solu- 
tion. On a small spatial scale the outflow expands and the 
Lorentz factor reaches a maximum of 7 = 2.8 in the inter- 
mediate region. This Lorentz factor is obviously smaller by 
a factor of 2 to 3 from observed values. We could rescale 
the whole solution to attain higher Lorentz factors but at 
the prize of high temperature in the jet. Similarly to the 
Solar wind this would require extra sources for the heat- 
ing to account for it. However we were mostly interested in 
describing the general features of the solution and try to 
explain the dichotomy between FRI and FRII jets rather 
than modeling specific observed jets. 

Then the solution obtained for a typical FRI jet recolli- 
mates, undergoes several oscillations and decelerates down 
to a Lorentz factor of 7 = 2 in the asymptotic region. This 
decrease is related to thermal compression of the jet in the 
asymptotic region by the outer medium. We insist on the 
fact that this recollimation occurs on a scale smaller than 
one parsec and that the external pressure is certainly the 
pressure of the surrounding disk wind rather than the exter- 
nal gas from the host galaxy. However if we assume that the 
extra pressure of the host galaxy can enhance the pressure 
in the disk wind, it is striking to note that our simple toy 
model for the spine jet shows a typical feature of FRI jets. 
Indeed FRI jets show a deceleration on the kiloparsec scale 
down to non relativistic speeds sometimes though they are 
usually highly relativistic on smaller scales. Besides that 
they are also known to have a rich ambient medium and 
that t he external gas pres sure is important as seen in the 
X-ray ( |Capetti et al |2002[ ). Finally the fact that the FRI jet 
radius is larger on the kiloparsec scale can be due precisely 
to this recollimation effect which enhances the density such 
that probably a larger part of the radio jet is emitting a s 
suggested in the double component jet of Sol et al. (19891. 



In the case of the FRII jets, the optimal value found is 
7/ = 0.86, which also gives a maximum Lorentz factor of 
7 ~ 3. The fact that r\ changes very little when we pass 
from the "FRI" -type solution to the "FRII" -type, shows 
that the outflow properties are highly sensitive to the rota- 
tion velocity of the corona. The solution obtained is how- 
ever very different as the gas expands monotonically and re- 
mains highly relativistic up to large distances, as observed 
in FRII jets. Moreover it is self collimated by its own mag- 
netic field throughout the length of the jet, something again 
that is characteristic of FRII jets wherein there is evidence 
that the host galaxy gas is rather underdense. 

On the other hand, the fact that in these two cases, 77 
is slightly smaller than 1 shows that the central launching 
region of the spine jet has to be slightly sub-Keplerian. 

In the present paper we did not discuss the case of 
Seyfert galaxies which are known to have outflows, though 
their winds are not very well collimated and are not rela- 
tivistic with velocities of the order of 30,000 km/s, as ob- 
served by the HST. Again, such flows can be understood in 
the frame of this simple self-similar model as non collimated 
solutions which are radial an d exist only if the vel ocity does 
not reach relativistic values ( Meliani et al.|200 6a). Such so- 
lutions are obtained if the magnetic rotator efficiency is 



very low, i.e. for e very negative, something similar to the 
solar wind. 

Altogether then, we may conclude that using a simple 
toy model for spine jets, the usual classification of radio 
sources can be understood on one hand by projection ef- 
fects and Doppler boosting and on the other by consider- 
ing the efficiency of the central magnetic rotator. Of course 
this does not exclude any of the other explanation such as 
the role of the external confinement or shear instabilities. 
Indeed our FRI-type solution is at least partially confined 
by the disk wind which may be a signature of external pres- 
sure confinement as well. This idea needs to be further ex- 
plored with more sophisticated models or simulations com- 
bining the central coronal jet with an external disk wind, 
something worth to pursue in another study. 

Next step under consideration is the extension of the 
model to the Kerr metric. This is a unique chance to con- 
struct the first analytical models for a jet around a rotat- 
ing black hole. One issue of this extension is to test the 
Blandford-Znajek mechanism (Blandford & Znajek 1977). 
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